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SUMMARY 


The application of the power cepstrum to measured helicopter-rotor acoustic data 
is investigated* A previously applied correction to the reconstructed spectrum, 
which requires prior knowledge of the echo amplitude or path length, is shown to be 
incorrect. For an exact echoed signal, the amplitude of the cepstrum echo spike at 
the echo delay time is linearly related to the relative amplitude of the echo in the 
time domain. It is shown that if the measured spectrum is not entirely from the 
source signal, the cepstrum will not yield the desired echo characteristics and 
cepstral aliasing may occur because of the effective sample rate in the frequency 
domain. In addition, the spectral analysis bandwidth must be less than one-half the 
echo ripple frequency or cepstral aliasing can occur. The power cepstrum editing 
technique is a useful tool for removing some of the contamination because of acoustic 
reflections from measured rotor acoustic spectra. Although the cepstrum editing 
yields an improved estimate of the free-field spectrum, the correction process is 
limited by the lack of accurate knowledge of the echo transfer function. An alter- 
nate procedure is proposed which allows the complete correction of a contaminated 
spectrum through the use of both the transfer function of the echo process and the 
echo delay time, but which does not require editing in the cepstral domain. 


1 . INTRODUCTION 

The concept of the cepstrum was first reported by Bogert, Healy, and Tukey 
(ref. 1 ) as a technique for detecting echoed signals in measured time series. Two 
general cepstral techniques have been used, the power cepstrum and the complex cep- 
strum. The power cepstrum is defined as the inverse Fourier transform of the natural 
logarithm of the power spectrum of the signal. The complex cepstrum is defined as 
the inverse Fourier transform of the complex logarithm of the Fourier transform of 
the signal. Both techniques make use of the convolution property of the Fourier 
transform, that a convolution in the time domain can be transformed to a product in 
the frequency domain. Both techniques also take advantage of a property of the loga- 
rithm function, that the logarithm of a product of two terms is equal to the sum of 
the logarithms of each term alone. One can consider a directly received signal and 
its echo as a signal convoluted with itself. The power spectrum of this time series 
is the product of the spectra of the direct signal and its echo. Taking the natural 
logarithm of the spectrum then separates the product into the sum of two natural log- 
arithm spectra, that of the direct signal and that representing the echo convolution. 
The unique characteristics of the individual cepstra of the direct signal and that of 
the echoes can allow editing of the echo contributions in the total cepstrum. Ide- 
ally, if the echo cepstrum can be removed, a transformation back to the frequency 
domain obtains the spectrum of the original direct signal. 

Although the power cepstrum technique allows the reconstruction of an uncontami- 
nated power spectrum, it does not maintain phase information, so that the original 
uncontaminated time series may not be reconstructed. To obtain the original time 
series, the complex cepstrum must be used. This technique is complicated by the use 
of the complex logarithm. Since the complex logarithm is a multivalued function, 
problems of phase discontinuities arise and the phase function must be "unwrapped" 
before the calculations can proceed. 


Over the last 20 years , the power cepstrum and complex cepstrum techniques have 
been applied to a wide variety of problems: seismic data (ref. 2), underwater sig- 

nals (ref. 3), voice analysis (refs. 4 and 5), brain wave signals (electroencephalo- 
graph (EEG) and visual evoked responses (VER) (refs* 6 to 8)), gearbox fault analysis 
(ref. 9), jet noise measurements (refs. 10 to 13), hydrophone calibration (ref. 14), 
reflector transfer functions and reflection coefficients (ref. 15), and internal- 
combustion-engine exhaust noise (ref. 16). References 17 and 18 give excellent 
reviews and bibliographies of the previous applications of the power and complex 
cepstra in addition to good reviews of the theoretical and practical considerations 
for each method* 

The objective of this work is to investigate the application of the power cep- 
strum technique to the analysis of model-scale helicopter-rotor noise spectra. 
Model-scale wind tunnel tests designed to study the acoustic radiation of helicopter 
rotors have been underway over the past two decades* Because of the relatively large 
size of model-scale rotors (e.g., a 1 /5-scale rotor typically has a 3-m diameter), a 
wind tunnel test must be conducted in a fairly large test chamber. In order to ac- 
quire high-quality acoustic measurements, the tests should ideally be performed in an 
anechoic facility, one in which the test chamber is treated to minimize acoustic re- 
flections which would contaminate the desired measurements. Unfortunately, there are 
few large-test-section wind tunnels in existence which are both anechoic and suitable 
for model rotor testing. For this reason, model tests are often performed in tunnels 
which have less than ideal acoustic qualities, and thus the measured data are contam- 
inated to some degree with acoustic reflections, in addition, acoustic measurements 
are often made outdoors at airports and heliports, and are thus easily contaminated 
with reflections. The power cepstrum technique has been successfully applied by two 
independent researchers to jet noise spectra which were contaminated with ground 
plane reflections (refs. 10 to 13). Although the technique was successful, this 
paper reports a minor error in the analyses employed in these previous efforts 
(refs. 10 and 12). However, the characteristics of jet noise and rotor noise are 
quite different, as there is a strong low-frequency tonal content in the rotor signal 
whereas jet noise is considered more broadband or random in nature. 

A joint experiment was performed in 1983 in the Langley 4- by 7-Meter Tunnel by 
NASA, U.S. Army, and Sikorsky Aircraft Company personnel (refs. 19 and 20). This 
experiment was the first phase of a three-phase program to investigate the interac- 
tion noise of main rotors and tail rotors. During the first phase, the noise of a 
main rotor alone was studied. The Langley 4- by 7-Meter Tunnel was not designed to 
be an anechoic test facility, and although measures have been taken to improve its 
characteristics, the acoustic data acquired contain contamination from the test cham- 
ber floor and walls. The acoustic qualities of the tunnel environment are addressed 
in references 20 and 21. The motivation for the research presented herein was to 
develop and validate a data reduction technique which could improve the accuracy of 
the measured acoustic data by removing the undesired reflections. 


2. THE POWER CEPSTRUM 


Theory 

One ideal reflection .- Recall the definition of the forward Fourier transform of j 
a general time function x(t) is 1 
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X(f) = F{x( t) } = I x(t) e' l27Tft dt 


and the inverse Fourier transform of X(f) is 


x(t) = F‘ 1 {x(f)} = I X(f) e l2lTft df 


To be an exact Fourier transform pair, the integration limits of equations (1) and 
(2) are infinite. When dealing with a finite amount of discrete data, however, the 
time integration must be performed over the period T and the frequency integral is 
valid only to the Nyquist frequency f N . 

The model of an acoustic signal g(t) with one ideal reflection is 


x(t) = g(t) + ag( t - t e ) 


where a is a constant and 0 < a < 1 .0. As discussed in section 5, the only con- 
straint on the nature of g(t) is that it have a continuous spectrum in the fre- 
quency range of interest. The Fourier transform of equation (3) is 


X(f) - ( [g(t) + ag( t - t^)] e i2lTft dt 

J -CO 


-i2TTft , I x -i2irft 

g(t) e dt + a I g(t - t^) e dt 


Using the substitutions 0 = t and 0 = t - t in each integral and the definition 


G(f) = I g( 0) e~ l27Tf0 d0 


we find that 


( -i2iTft \ 

= ll + a e e /< 


The power spectrum of x(t) may be defined as 
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X(f ) I 2 = X(f )X*(f ) 



-i2irft 

+ a e 



+ a e 


G(f )G*(f ) 


1 + 


2a 


1 + a 


— cos(2nft ) (1 + a ) 
2 e 


G(f) | 2 


( 6 ) 


Note that in the absence of a reflection (i.e., for a = 0), |x(f)| = | G ( f ) | • 

Thus, the terms of the first two brackets of equation (6) are due to the presence of 
the reflection (i.e., the echo). The cosine term causes the characteristic ripple 
seen in the power spectra of data containing a reflection. To decouple the terms of 
equation (6), the natural logarithm is used: 


In X(f)r = In 


1 + — — cos(2TTftJ + ln(1 + a 2 ) + In | G ( f ) | 2 
1 + a 


(7) 


To further decouple equation (7), the power series expansion of the natural logarithm 


In ( 1 + x ) 




which is valid for |x| < 1, is applied to the first term of equation (7) because, 
for 0 < a < 1.0, 


2a 


1 + a! 


- cos(2TTft ) 
2 e 


< 1 


Further simplification of the series expansion is obtained by letting 


4> = 


2a 


1 + a 


Thus, equation (7) is written as 


ln|x(f)|^ = ln|G(f)j 2 + ln(1 + a 2 ) + £<j> cos(2irft e ) 

- 4- 4> 2 cos 2 ( 2 tt f t ) + 1 <t> 3 COS 3 (2TTft ) + ...1 
2 e 3 e J 


( 8 ) 
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The power cepstrum is defined as the inverse finite Fourier transform of the 
natural logarithm of the power spectrum: 


C (x) = F _1 {in | X(f ) I 2 } 
P 


(9) 


Note the similarity to the autocorrelation function 

v t) = i x(f) i 2 } 


Through substitution of equation (8) into equation (9) and use of identities relating 
(cos 0 ) n to cos(n6), the power cepstrum is found to be 




C (t ) = I N ln[ | G ( f ) | 2 ] 


N 


2 1 i2nf X ,, I “N , 2. i2irfx 

e df + | ln(1 + a ) e df 


4> f N cos(2irft ) e ^ 2ir ^ T , 3 f _ 1. (j, 2 f N l-[ i + cos(4irft )] e 

J - f e 2 J-f 2 e 

N N 


i2irfx 


df 


1 ,3 f f N 1 

J. f * 

N 


[3 COS(2TTft ) + cos(6irft )] e ^- 2lT ^ T df + 
L e e J 


( 10 ) 


A lengthier expansion of this equation, to the tenth power of (J) cos(27rft e ), is pre- 
sented in appendix A for completeness and clarity. Here, the first term is simply 
equal to the power cepstrum of the direct signal g(t). To perform the 

py 

integration of the remaining terms, note that each can be reduced to the integral of 
cos(2Trnft e ) e : *' 27TfT (where n = 0, 1, 2, ...) and recall that 



cos(n2Trft ) 
e 


i2irf T 
e 


df 


Jsin[(nt e - T)2irf N ] sin[(nt e + t)2irf N ] 


(nt - x)2irf 
e N 


(nt + x)2irf 
e N 


(n=0, 1, 2, ...) 


It is interesting to note that if the integration could be performed over infinite 
limits, the (sin x)/x functions would become delta functions. 
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Now equation (10) may be simplified to 


C (x) = C (x) + 2f [ In ( 1 + a 2 ) - a 1 
p pg n l 0 j 


sin( 2irf ,t)' 
N 


2*V 


00 ( sin [ ( nt - x ) 2irf 1 sin[(nt + x)27rf 1 

- y \ L e NJ_ L e N_ J 

+ N " a n l (nt - t ) 2irf + (nt + x)2iTf 

n=1 l e N e N 


(11 ) 


The second term represents a peak at x = 0, and the third term represents an infi- 
nite series of peaks (called "rahmonics" ; analogous to harmonics of a frequency) at 
multiples of the echo delay time t 0 with amplitudes (f N a n ) which are related to the 
ratio a. 


The typical shape of the power cepstrum of the direct wave (C (x) for a gen- 
eral random signal) is shown in figure 1 . It is similar in shape to the autocorre- 
lation of a random signal, as we would expect since the cepstrum would be identical 
to the autocorrelation if the natural logarithm were not used. The presence of 
reflections in the signal causes the series of spikes (the (sin x)/x series) in 
the cepstrum. To remove the reflections, one can remove the spikes at ±nt e 
(n=1, 2, 3, ...) by interpolating new values based on the adjacent values. 

2 

It can be shown that the coefficient ln(1 + a ) - a Q of the second term in 
equation (11) (the spike at x = 0) is identically zero. (See appendix A.) Thus, 

the only contribution at x = 0 is from C_( x ) . it should be noted that Syed et al. 

Py 2 

(ref. 10) used a different series expansion for ln(1 + a ) and approximated the 
series to only the fourth power. He found that the spike at x = 0 was multiplied 
by ln(1 + a 2 ) - 4.35 [ ( 4> 2 /4 ) + ( 3<J> 4 /32 ) ] and that this additional term should be 
removed. However, since there is a large contribution at x = 0 from C^g ( X ) , it is 
difficult to edit the cepstrum at that location. Thus, he concluded that one should 
not edit the cepstrum at x = 0; rather, one should edit only the spikes at 
i = ±nt e , and when tlie cepstrum is transformed back to a power spectrum, the addi- 
tional x = 0 term because of the echo becomes a constant amplitude correction to be 
applied to the reconstructed spectrum. This correction would require knowledge of 
a. However, since this term is actually zero, no amplitude correction is required. 


In addition, work reported by Miles, Stevens, and Leininger (ref. 12) concludes 
that editing of the cepstrum spikes results in "a spectrum due to intensity addition 
only of the direct and reflected signals.” They also found that the coefficient of 
the x = 0 spike is related to ln(1 + a ) and that knowledge of a would allow a 
correction to obtain the free-field spectrum. This also is an incorrect conclusion. 

Since the second term in equation (11) is actually multiplied by zero, the 
actual cepstrum is 


C (x) 
P 


C ( 

pg 


00 ( sinf (nt - x)2TTf 1 sin[(nt + x)27Tf 1 

t n + f 1 L e L e N__ J 

N “ a n| (nt - x)2*rrf + (nt + x)2irf 

n=1 ^ e N e N 


( 12 ) 
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After the ±nt e spikes are removed from the cepstrum, the new edited cepstrum 
C^(x) is 


C' (x) = C (x) = F -1 {ln|G(f) | 2 } 
p pg 1 


To find the free-field power spectrum, a Fourier transform is again performed 


In | G( f ) | 2 * F{c^(x) } 


and the exponential function is used 


|G(f ) | 2 = e 


2 _ _ F { C p< T 0 


(13) 


This results in reconstruction of the free-field power spectrum of the original sig- 
nal and does not require any amplitude corrections or prior knowledge of the echo 
amplitude. 


Though the cepstrum spikes yield information on the delay time of echoes, they 
may also be used to estimate the relative amplitudes of the echoes. Through the use 
of equation (10), it can be shown that the coef f icient^of ^ the first echo spike of the 
cepstrum (i.e., the coefficient of the / cos(2itfx e ) e df terms) is f N a 1 , 

where 


1 3 ,3 1 10 ,5 1 35 J 1 126 x 9 


a i = * + 34* + S 76 * + 7 64 * + 9 


256 


♦ 


(14) 


with 


♦ = 


2a 


1 + a 


It can be shown that this is a convergent series for 0 < a < 1 .0 (see appendix B) , 
and the series can be rewritten as follows: 


a 


1 


» 2i —1 

V 9 3 ( 2j - 2)1 

2i -2 

j-1 2 Z1 ^(j - 1 ) I j ! 


2a 


(15) 


If the value of C (x) is small at t e (the value of x for which the first echo 

spike occurs, i.e., x = t g ), then the amplitude at t e should be essentially 2f N a. 

For many functions, C (x) typically decreases very quickly with increasing x, so 

py 
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often this amplitude is valid. Thus, given the cepstrum of a signal, the amplitude 
of the first echo spike may be used to find the value of the relative amplitude of 
the reflection. 

Two ideal reflections .- The model of an acoustic signal with two ideal reflec- 
tions is 


x(t) = g(t) + ag(t - t 1 ) + 3g(t - t 2 ) 


(16) 


Using the same approach as for the one-echo case (refer to eg. (6)), we may express 
the power spectrum of x(t) as 


|x(f ) | 2 


2a 


(1 + 


a 2 + 


6 2 ) 


cos ( 2 Trf t ^ ) + 


23 


(1 + 


2 

a + 


e 2 ) 


cos( 2nft^) 


+ — cos[2Trf(t - t )]l(1 + a 2 + 0 2 ) |G(f ) | 2 

(1 + a + 0 ) J 


(17) 


To decouple the terms of equation (17), the logarithm is used and is approximated by 
the power series expansion to the third power to produce 


ln|x(f ) | 2 - In | G ( f ) | 2 + 2a cos(2irft ) + 20 cos(2Tfft 2 ) 

- a 2 cos(2irf2t 1 ) - 0 2 cos ( 2irf 2t 2 ) + j a 3 cos ( 2trf 3t 1 ) 

+ j 0 3 cos(2irf3t 2 ) - 2a0 cos[2irf(t + tj)] 

+ 2a0 2 cos[2irf(t i + 2t 2 ) ] + 2a 2 0 cos[2irf(2t 1 + t 2 ) ] 
This approximation is only valid when 
| (a + 0) 2 + 2(a + 0) | <1 
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Performing the inverse Fourier transform to obtain the power cepstrum results in 


C (T ) * C (T) + f 
p pg N 


n=1 


n 


"sin[(nt 1 ± T)2irf N ] 

(nt ± t )2irf 
1 N 


* |'sin[(nt ± x)2irf N n • /sin{[n(t + t ) ± T]2irf N }' 

+ 2- b - S TZZ -r-TT^ rz 1 + 2- c 


. n I (nt ± t ) 2ir f 
n=1 L 2 N 


n=1 n \ [n(t 1 + t 2 ) ± T]2irf N 


J d / Sin ^ n(t 1 + 2t 2 } * T ^ 2ltf N^ £ ^sin{[n(2t 1 + t 2 ) ± T]2*f N }' 


n=1 n \ + 2t 2 ) ± t]2irf N / n =1 n \ [n(2t., + t 2 ) ± T]2irf N 


(18) 

2 2 

The factor 1 + a +3 in equation (17) is cancelled by the aQ term from the 
power series expansion of the logarithm, as is found for the one-echo case (appen- 
dix A). No attempt has been made to deduce the exact converged values of the coeffi- 
cients a n , b n , c n , ... (n > 0) as is done for a-j of the one-echo case. However, 
this power series approximation shows that 

2 2 3 2 

a 1 « 2a, a 2 * -a , a « -j a , b * 23, b * -3 , 

b 3 » j 3 3 f c * -2a3# d » 2ot3 2 , and e 1 » 2a 2 3 


For the two-echo case, the dominant echo spikes occur at t-j , 1 2 , 2t<, , 2t 2 , and 

t 1 + t 2 r and the amplitudes at t 1 and t 2 are the same as for a one-echo case. No 
spikes occur at differences in delay times. The rahmonic spikes (those at sums of 
multiples of t 1 and t 2 ) are of lesser amplitude, since they are multiplied by pow- 
ers of a and 3 and in practical cases a and 3 are less than unity. Thus, in 
analogy to the one-echo case, the cepstrum echo-spike amplitudes are linearly related 
to the actual echo amplitudes. 


Numerical Examples 

Direct signal .- To illustrate the results of the preceding analysis, a numerical 
simulation of a transient signal was performed. The sample transient signal has the 
exponential form 


x(t) 


-k t-t 
e 1 


(19) 


The time history, power spectrum, and power cepstrum are shown in figures 2(a) to 
2(c). The power cepstra presented in this report are normalized by twice the Nyquist 
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frequency (2f N ), so that the amplitude of the first cepstrum spike represents the 
relative echo amplitude as shown by equation (15). In addition, the first few points 
near t = 0 are not plotted to increase the plot resolution for the cepstrum at 
large values of T. 

One ideal reflection.- The sample signal with one echo has the form 


x( t) 


-k 

e 



-k 

+ a e 



( 20 ) 


This signal and its power spectrum for a - 0.5 and = 0.0018 sec are shown in 
figures 3(a) and 3(b). Note the large-amplitude ripple induced in the spectrum, and 
that the overall level is 1 dB higher than the original signal (fig. 2(b)). Fig- 
ure 3(c) shows the normalized cepstrum of the signal of equation (20). The quickly 

decaying C (t) function is shown near x = 0, and the spikes at multiples of the 
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echo delay time are clearly shown. The amplitude of the first cepstrum spike of 
figure 3(c) is 0.50 (with the neighboring cepstrum values nearly zero), which illus- 
trates the equality with the relative echo amplitude of 0.50 proven by equation (15). 
These spikes are easily removed by linearly interpolating new values based on adja- 
cent cepstrum values. In this case the editing procedure is fairly straightforward, 
since the cepstrum levels are very small for x > 0.001 sec, and the data are 
numerically "clean • " 

The reconstructed power spectrum is shown in figure 3(d) and is simply the expo- 
nential of the forward Fourier transform of the edited cepstrum. This reconstruction 
is identical to the Fourier transform of the isolated wavelet (fig. 2(b)) and it is 
clear that the previously used amplitude corrections (refs. 10 and 12) need not be 
applied to this result. The difference between the contaminated spectrum (fig. 3(b)) 
and the edited spectrum (fig. 3(d)) is shown in figure 3(e), which clearly shows the 
sinusoidal nature of the echo contamination. 

Two ideal reflections .- To demonstrate the application of the technique to a 
signal with two echoes, figure 4(a) shows the same wavelet as shown in figure 2(a) 
with two echoes. Figures 4(b), 4(c), and 4(d) show the power spectrum, power cep- 
strum, and reconstructed "f ree-f ield" power spectrum of this signal. As shown in 
equation (18), the effect of multiple echoes is to create cepstrum spikes not only at 
multiples of the echo delay times but also at sums of multiples of the echo delay 
times. This effect is shown in figure 4(c). The amplitudes of the cepstrum spikes 
at x = 0.0018 and 0.0068 sec correlate well with the expected values of 

fa « f 2a and f b » f 23 found from equation (18). The amplitudes of the cep- 
N 1 N N 1 N 

strum spikes at multiples and sums of multiples of the two echo delay times also 
correlate well with the approximated values of their coefficients. (Recall that the 
cepstrum is normalized by twice the Nyquist frequency.) Again, since this is a 
numerically "clean" example, editing of this signal is fairly straightforward. 

Twenty- two of the highest level spikes in the cepstrum of figure 4(c) were edited 
out. Figure 4(d) shows the resulting edited f ree-f ield power spectrum, a good 
reconstruction of the free-field spectrum (fig. 2(b)). 
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3. EFFECT OF ADDITIVE NOISE 


Theory 

An important case occurs when noise is present in addition to the direct and 
reflected signals. The model for this is 


x ( t ) = g(t) + ag(t - t ) + n(t) 


The Fourier transform is 


-i27rft. 


G ( f ) + N(f) 


The power spectrum may be expressed as 


X(f ) = X(f )X*(f ) 


[ 1 + — 2a cos(27rft )1 (1 + a 

1 + a 2 e J 


2 ) I G( f ) I 2 


+ | N( f ) | 2 + 2 { [l + a COs(2irft e ) ] [G f (f )N r (f ) + G^(f)N^(f)] 

- a sin(2irft ) [g (f )N. (f) - G. (f)N ( f ) ] } 
e L r i x r J J 


where G r , N r , G^, and N.^ are the real and imaginary parts of G(f) and N(f). 
The first term is identical to that of equation (6), and the last three terms are due 
to the noise. If N(f) is small compared with G(f), the contribution of the noise 
terms may be negligible, but if N(f) is large enough, these terms will contaminate 
the power spectrum and^the power cepstrum. The contamination consists of the noise 
power spectrum | N ( f ) | and the "cross terms" (i.e., G r (f)N r (f), G^(f)N^(f), 

G r (f)N i (f), and G i (f)N r (f)). 

A probabilistic study (refs. 22 and 23) has concluded that the effect of the 
cross terms is to lower and broaden the cepstrum echo peaks. The degree of con- 
tamination is related to both the frequency dependence of the signal-to-noise ratio 
and the relative bandwidths of the signal and noise (refs. 22 to 24). There are 
several approaches which can be used to avoid or decrease this noise contamination. 

If n(t) could be measured without the source signal present, it could be subtracted 
from x(t). Similarly, if N(f) is known, it could be removed from X(f) or 
| X ( f ) | 2 • Ensemble averaging in the time domain would also serve to remove some of 
the noise content, but such time averaging would also remove any random content of 
x ( t ) , which may or may not be desirable. Averaging of the Fourier transforms would 
yield essentially the same results as averaging in the time domain. Averaging of the 
power spectra would eliminate the noise cross terms, since the signal and noise func- 
tions are incoherent. However, averaging of the power spectra will not remove the 
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noise power spectrum | N ( f ) | 2 • Averaging in the cepstral domain has also been sug- 
gested (ref. 15). Such averaging procedures and data manipulation must be carefully 
chosen based on the character of the data (i.e., random, coherent, or both). 


Numerical Examples 

To qualitatively investigate the effect of additive noise, a composite signal 
consisting of the exponential impulse and echo of figure 3 with an additive random 
"background” noise function 


x ( t ) 


-k t-t 


-k 

+ a e 



+ n(t) 


was constructed to simulate the power spectrum described by equation (21). (See 
fig. 5.) The noise function n(t) was obtained from a software random number gener- 
ator. The signal-to-noise ratio of the exponential signal and echo to the noise is 
about 20 dB. The amplitude of the cepstrum spike at x = 0.0018 sec is signifi- 
cantly less than expected for an echo of ot = 0.5, and rahmonics are not present be- 
cause of the noise, as shown in figure 5(c). The cepstrum spike amplitude is approx- 
imately 0.10, as opposed to an amplitude of 0.50 shown in figure 3(c) for the "clean" 
transient signal. Removal of this single spike results in a partial improvement in 
the echo ripple in the power spectrum. (Compare fig. 5(d) with the corrected spec- 
trum without noise, fig. 3(d).) This illustrates that the additional noise terms of 
equation (21) can contaminate the cepstrum sufficiently to destroy the effectiveness 
of the approach, even if the noise N(f) is small compared with the desired signal 
and echo. This example can be considered a "worst case" for two reasons: (1) no 

noise-removal approach, such as temporal, spectral, or cepstral averaging, is at- 
tempted; and (2) the relative signal-to-noise ratio of the exponential signal to the 
noise is quite small over most of the spectrum above a frequency of 10 kHz. 

To verify the results of references 22 to 24 (i.e., that the noise contamination 
is related to the frequency variation of the signal-to-noise ratio and to the rela- 
tive bandwidths of the signal and noise), two other cases were investigated. The 
first was the case of increasing the signal-to-noise ratio over that of figure 5; the 
second was the case of increasing the signal-to-noise ratio and widening the signal 
bandwidth relative to the noise bandwidth. 

Figure 6 illustrates the effect of increasing the signal-to-noise ratio of the 
previous example from 20 to 32 dB. The cepstrum shown in figure 6(c) is very similar 
to that of the previous case, except that the spike amplitude has increased by a fac- 
tor of two to a value of 0.22. This is closer to the desired value of 0.50 which 
would result if no noise were present. If we edit the cepstrum and perform the in- 
verse Fourier transform, the corrected spectrum of figure 6(d) is found. This is a 
better estimate of the free-field spectrum than figure 5(d), although a small echo 
ripple remains. 

Figure 7 illustrates the case of a larger signal bandwidth relative to the noise 
bandwidth. in this example the exponential factor k was increased from 6000 to 
20 000, providing a broader signal spectrum. The normalized cepstrum (fig. 7(c)) 
exhibits a prominent echo spike at the correct delay time, with an amplitude ap- 
proaching the desired level of 0.50. The first echo spike and one rahmonic were 
edited and the reconstructed spectrum is shown in figure 7(d). 
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4. EFFECT OF ECHO DISTORTION 


Theory 

The case of a direct signal and a distorted echo is also of great practical 
importance and can be modeled as 


x ( t ) = g(t) + h(t - t e )*g(t) 


where the asterisk denotes a convolution and h(t) is the impulse response function 
of the echo process. The Fourier transform and power spectrum may be written as 


X(f ) 


[l + H(f) e 


- 2 * ft e] 


G(f) 


|x(f) 


= Li 


+ H(f) e 
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( 22 ) 


If we take the logarithm and use a power series expansion, the log power spectrum is 
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where H(f) is the transfer function of the echo process and the superscript aster- 
isk denotes the complex conjugate. Note that if H(f) is a constant, causing a 

sin[(nt e - x)2irf N ] 

perfect echo, then h(t - t„) reduces to the — -r— r— — term of 

, i / h <i \ ( nt ™ T ) 2 TT t 

equation (11). e N 


This case has been studied in detail by Hassab and Boucher (refs. 22 and 23) and 
by Bolton and Gold (ref. 15). The conclusions of reference 22 show that with distor- 
tion, the cepstrum echo spikes are lowered and broadened, qualitatively a similar 
effect as was shown in the presence of additive noise. This effect on the cepstrum 
spikes depends on the transfer function H(f) of the echo process, which of course 
is generally unknown. The echo process can be considered as the sum of the reflec- 
tor’s transfer function and the acoustic propagation path. 

The echo information in the cepstral domain is now "windowed," or smeared, by 
the transfer function in the same manner that windowing in the frequency domain 
occurs. The echo no longer appears as a discrete spike at the echo delay time and 
editing the cepstrum to remove the echo becomes substantially less straightforward. 


Correction of Echo Transfer Function 

To remove the contamination in the power spectrum due to the distorted echo, 
another approach is suggested. This approach requires an accurate knowledge of the 
transfer function of the echo process and an accurate knowledge of the echo delay 
time t e . The complex impedance of the reflector and the acoustic propagation ef- 
fects can be obtained through the use of standard techniques. The complex impedance 
can also be found with the cepstrum technique of reference 15. An accurate estimate 
of the echo delay time can be obtained easily from the power cepstrum of the contami- 
nated signal. Now equation (22) may be rewritten as 


| X( f ) I 2 = [l + 2H r (f) cos(2irft ) + 2H i (f) sin( 2irf t g ) + | H( f ) | 2 ] | G( f ) | 2 (23) 

Rearranging these terms, we can find the desired free-field power spectrum | G ( f ) | 2 
by correcting the contaminated power spectrum |x(f)| as follows: 


Mf )| 2 ^ 5 

1 + 2H (f) cos(2irft ) + 2H. (f ) sin(2irft ) + | H ( f ) | 
r e i e ' 


The advantages of this technique are that the case-by-case editing of the cepstrum 
and the inverse transform processes are not required and that a routine correction to 
all measured power spectra can be employed. Elimination of the cepstrum editing and 
the inverse transform processes will be particularly beneficial when the cepstrum is 
difficult to edit due to noise or periodic content. The disadvantage to this tech- 
nique is that it requires a knowledge of the echo transfer function H(f). 
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Example of Transfer-Function Correction 

To verify the validity of this transfer-function correction procedure , the pro- 
cedure was applied to measured acoustic impulse data containing an echo. The impulse 
data were obtained during a model helicopter-rotor wind tunnel test in the Langley 4- 
by 7-Meter Tunnel in which remotely ignited blasting caps were mounted in the center 
of the tunnel test section. The direct blast signal and echoes were measured at all 
the microphone locations in the test section. Since the direct and echoed signals of 
the impulse occupy different regions of the acoustic time series f the Fourier trans- 
form of each impulse (X^(f) and X e (f)) may be calculated separately, and then the 
echo transfer function may be estimated. 

Equation (22) may be expressed as the sum of the Fourier transforms of the 
direct signal (X d (f)) and the echoed signal (X e (f)) as follows: 

[ -i2iTft 1 

X(f) = Li + H(f) e e jG(f) = x d (f> + X e (f) 


where 


X d (f) = G(f) 


and 


-i2irft 


X e (f) * H(f) e 


G ( f ) 


The ratio of X d (f) and X e (f) are used to define a function R(f): 


R(f ) = 


X (f) 
e 

x,(f ) 

d 


(25) 


With R(f), H(f) may be found from 


i2irft 

H(f) = e 6 R(f) 


From equation (24) with H(f) and t Q , the free-field spectrum may be found. 

Figure 8(a) presents the impulsive acoustic signal with a relatively strong re- 
flection (echo) measured with a microphone very close to the tunnel test section 
floor. The regions marked in the figure indicate the data which were used to inde- 
pendently calculate the Fourier transforms of the direct and echoed signals. The 
same record length was used for each, with all nonpertinent data points before or 
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after the transient impulses set to zero. Figures 8(b) and 8(c) show the power spec- 
tra of the direct and echoed signals. These data were then used to calculate R(f) 
from equation (25). Figure 8(d) shows the power spectrum of the direct plus the ech- 
oed signal. The objective is to correct this contaminated spectrum for the nonlinear 
echo and to reconstruct the uncontaminated spectrum of the direct wave (fig. 8(b)). 

To calculate the function H(f), an estimate of the echo time delay t 0 is re- 
quired. This value is estimated to be 0.003 sec from the power cepstrum of the con- 
taminated spectrum (shown in fig. 9(a)). The transfer^function H(f) may now be 
found from R(f) and t 0 . The power spectrum | H( f ) | is shown in figure 9(b). 

When the contaminated spectrum is corrected for H(f) as defined by equation (24), 
the spectrum shown in figure 9(c) is found. This spectrum compares very well with 
that of the direct wave presented in figure 8(b). 

It is interesting to compare this correction procedure with the process of 
editing the cepstrum. The cepstrum shown in figure 9(a) was edited by setting all 
values greater than t = 0.002 sec to zero. The resulting spectrum is shown in 
figure 9(d). (Such severe cepstrum editing is valid when the time signal is tran- 
sient or nonperiodic, since the cepstrum of the direct signal is small for large 
values of t. The hump in the cepstrum near t = 0.003 sec (fig. 9(a)) is attrib- 
uted to the windowing effect of the nonlinear transfer function acting on the echoed 
signal.) In a comparison of the two reconstructed spectra with the spectrum of the 
direct signal alone (fig. 8(b)), the spectrum corrected for the echo transfer func- 
tion (fig. 9(c)) is clearly superior. The spectrum corrected with the cepstrum cor- 
rection process (fig. 9(d)) yields a satisfactory result only at frequencies less 
than 2 kHz. 

While performing these calculations, we noticed that the transfer functions cal- 
culated were extremely sensitive to the specific data used in each transform, the 
sample length, the sample rate, the transform size, and other particulars of the 
analysis. Although the spectral shapes were generally the same, the amplitudes would 
vary widely with small differences in the analysis approach. We feel that this was 
due to the paucity of data available to estimate the transfer function. Other 
approaches to finding the echo transfer function may prove more successful than that 
employed in this work. 


5. APPLICATION TO PERIODIC FUNCTIONS 

Much previous work on the topic of cepstrum echo removal has concluded that the 
technique is most successful if applied to cepstra which decay rapidly with quefrency 
t. Typically these cepstra are for random or noncoherent functions. Deterministic 
or periodic functions, which are coherent, generally exhibit periodic characteristics 
in their cepstra. (Recall that the autocorrelation of a sinusoid is sinusoidal.) In 
preparation for the application of the cepstrum echo removal technique to measured 
helicopter-rotor acoustic data (which is considered highly tonal noise), the applica- 
tion of the technique to a generalized periodic function was investigated. The exam- 
ple periodic function considered is a set of harmonic sine waves of equal amplitude. 
Several examples are considered which have different contributions (i.e., source- 
related, non-source-related, and leakage-related contributions) at the off-harmonic 
frequencies • 
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Numerical Examples 


Direct signal .- Figure 10(a) presents the power spectrum of the example periodic 
function , a series of 25 equal-amplitude harmonics of a fundamental 5-Hz tone with no 
echo contamination. The source signal has a constant off-harmonic level of about 
76 dB. Figure 10(b) presents the normalized power cepstrum of this signal. Note the 
strong periodic character of the cepstrum. The spikes at t = 0.20 (t f ) and 
T = 0.40 are due to the fundamental tone frequency of 5 Hz. The oscillatory period- 
icity of the cepstrum is due to the harmonics present. 

One ideal reflection (t 0 = 0.07 sec).- Figure 11(a) presents the power spectrum 

of the same signal as used for figure 10 with the addition of an echo (a = 0.5) at a 
delay time of 0.07 sec. The cepstrum (fig. 11(b)) clearly identifies this echo at 
x = 0.07 sec (t e ), a nd one rahmonic is evident at x = 0.14 sec (2t e ). Note that 
the magnitude of the echo spike at x = 0.07 sec is very close to the expected level 
of 0.50 (a = 0.5). Editing and interpolating of this type of cepstrum is not 
straightforward because of its harmonic character. Figure 11(c) presents the edited 
power cepstrum, and figure 11(d) presents the reconstructed power spectrum. A 
substantial part (though not all) of the echo ripple has been removed. 

One ideal reflection (t 0 = 0.07 sec) with leakage.- Another case of interest is 
that of the series of harmonic tones with no off-harmonic spectral content except 
that caused by numerical leakage or by the spectral leakage associated with use of a 
sample period which is not exactly an integer multiple of the fundamental period. 
Figure 12(a) presents the case of the same series of tones and an echo (a = 0.5 and 
t = 0.07 sec) as used for figure 11(a) with no other source-related energy in the 
off-harmonic bands. The leakage energy is apparent in between the harmonic tones at 
a level of about 10 dB. The power cepstrum (fig. 12(b)) contains an echo spike at 
the correct delay time, although the spike amplitude is somewhat lower (0.40) than 
the expected value of 0.50. This lowering of the echo-spike level is attributed to 
the lack of significant energy in the off-harmonic bands to carry the echo ripple 
information. 

One ideal reflection (t 0 = 0.29 sec).- Figure 13 shows the same periodic 

function presented in figure 12 with an echo of a = 0.5 at t e = 0.29 sec. 

The power cepstrum (fig. 13(b)) of this signal exhibits an echo spike at 
x = 0.09 sec. In this case the echo ripple in the power spectrum has effectively 
, been undersampled, and the information has been aliased in very much the same manner 
as aliasing in the frequency domain can occur. Instead of identifying an echo at 
x = t e (i.e., at 0.29 sec), the process has identified an echo at x = t 0 - t f 
(0.29 - 0.20 = 0.09 sec), an aliased time delay. (Recall that t f is the period of 
, the 5 -Hz tone.) This misidentif ication is attributed to the following: because 

there is not significant source-related energy in the off-harmonic bands to represent 
the echo ripple in the spectrum, the harmonic tone frequency becomes the "effective 
I sample rate" of the spectrum; and because the harmonic tone frequency (sample rate) 

| of 5 Hz is greater than the echo ripple frequency ( 1 / ti e » 3.5 Hz), aliasing occurs. 

This phenomenon is analogous to the sampling theorem for time series, which states 

: that to resolve frequencies with a period of T, one must use a sample interval of 

i less than T/2. Thus, to resolve a frequency of 3.5 Hz in the spectrum, the 

j effective sample rate should be less than 1.7 Hz. 

I One ideal reflection (t e = 0.07 sec) with additive noise.- A final example 

I (figs. 14 and 15) involves the same periodic function as that used for figures 11 to 
13 with an off-harmonic level which is unrelated to the source signal; this results 
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in an additive ambient noise level. This example, with an echo at t 0 = 0.07 sec, is 
presented in figure 14(a). The cepstrum (fig. 14(b)) only weakly exhibits the echo 
spike. The same source signal with an echo at t Q = 0.29 sec is presented in fig- 
ure 15(a). Again the cepstrum echo spike is very weak compared with the previous 
cases, and the echo time delay has been aliased from its actual value. 


Discussion 

The general conclusion from these results is that the measured source spectrum 
must be a continuous spectrum, that is, the spectral energy must all be due to the 
source signal and its echoes. If there are regions of the measured spectrum which 
are due to an extraneous source, then the echo ripple information will not be repre- 
sented at these bands and the detection of the echo will be severely deteriorated. 
Although the numerical examples may seem purely academic, the situations modeled can 
easily occur in experimental measurements. Consider, for example, the measurement of 
a source with an echo when the source has a highly tonal or periodic content. The 
source may also have broadband content at the nontonal frequencies. If the off- 
harmonic source energy is significantly lower than the tone energy, it may be hidden 
by instrumentation noise during the data acquisition process. Thus, one could be 
faced with data very similar to that presented in figure 14, a periodic direct signal 
with an echo at harmonic frequencies, and a non-source-related off-harmonic "noise 
floor" which would not represent the echo ripple. 

Other conclusions related to the application of the cepstrum to periodic data 
are as follows. The editing and interpolating procedure can be quite difficult, with 
potentially erratic results. For such signals, the correction technique presented in 
section 4 may be more useful. Clearly, the dynamic ranges of the data and of the 
data acquisition process are of importance, as is the origin of the spectral content 
(source- or non-source-related contribution). In addition, the ratio of the analysis 
bandwidth (1/T) to the echo ripple frequency (1/t e ) is important, in analogy to 
Shannon's sampling theorem, it is concluded that the analysis bandwidth should be 
less than one-half of the echo ripple frequency 1 /t e , to prevent aliasing of the 
echo time delay in the cepstral domain. This is a similar conclusion to that pre- 
viously reported (ref. 15), that the echo delay time should be less than one-half of 
the sample record length. 


6. APPLICATION TO HELICOPTER-ROTOR ACOUSTIC DATA 

The cepstrum correction technique was applied to measured model helicopter-rotor 
acoustic data containing reflections from a wind tunnel wall. These data were ac- 
quired in the Langley 4- by 7-Meter Tunnel. Acoustic results are reported in refer- 
ences 20 and 21 • A diagnostic test was also performed to assess the magnitude of 
possible reflections from an impulsive point source at each microphone and to iden- 
tify an approximate echo arrival time at each microphone. Estimates of the relative 
amplitudes were made from the acoustic time histories of these events, and estimates 
of the delay times were made from the power cepstra. 


Helicopter-Rotor Noise Model 

It should be recognized that helicopter-rotor noise consists of two general 
contributions. The highest noise levels are generated at a series of low-frequency 
harmonics of the blade-passage frequency f (equal to the number of rotor blades 
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times the shaft rotation frequency; typically, f Bp = 15 to 20 Hz). This acoustic 
signal is highly periodic and discrete. The other major source is helicopter-rotor 
broadband noise (ref. 25), and although this source has some energy at the lower fre- 
quencies (those less than 1 kHz), most of the energy is in the middle frequency range 
(i.e., 1 to 5 kHz). This contribution is a more random and broadband noise source. 
Thus, the helicopter-rotor noise signal can be modeled as follows: 


N 

x ( t ) = V a(nf ) sin( 2nnf t) + r(t) + n(t) 
BP BP 

n=1 


where a(nf 0p ) is the amplitude envelope of the helicopter-rotor harmonics, N is 
the number of rotor harmonics, ' r(t) is the rotor broadband noise signal, and n(t) 
is the additive wind tunnel background noise. The rotor signal with one perfect 
reflection can be modeled as 


x(t) = T a (nf ){sin(2irnf t) + a sin[2Trnf (t - t )]} 
BP 1 BP L BP e J J 

n=1 


+ r(t) + ar(t - t ) + n(t) 
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Application 

Power-averaged spectra .- The cepstrum technique was applied to the rotor 
acoustic data measured at a microphone located in the tunnel flow for a typical 
rotor test condition (i.e., 60 knots, 2° rotor angle of attack, microphone 3, and 
f Bp = 100 Hz). Since this was a 1/5-scale rotor model, the rotational speed was 
increased by five, causing the f B ^ value of 100 Hz. Fifty power spectra of the 
rotor signal were averaged to obtain a power-averaged spectrum as follows: 



This is the conventional method of spectral averaging, as it preserves all the acous- 
tic energy in the signal and does not remove random content. Another result of this 
spectral averaging process is to eliminate some of the cross term contributions of 
equation (21) as discussed in section 3. Note the ripple in the spectrum 
(fig. 16(a)), particularly obvious between 2000 and 6000 Hz. The period of this rip- 
ple is approximately 500 Hz, which corresponds to an echo time delay of 0.0020 sec. 
From the impulse data results, which are presented in table III of reference 20, we 
anticipate reflections occurring at echo delay times of approximately 0.0025 and 
0.0240 sec with relative amplitudes of 0.49 and 0.10. The cepstrum is shown in fig- 
ure 16(b), with two prominent spikes near 0.0020 and 0.0200 sec. (In these plots the 
first few points near zero have been set to zero to increase the plot resolution.) 

The rahmonics of these two spikes are not obvious • Note that except for the echo 
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spikes, the cepstrum shows a somewhat repeatable pattern. This pattern repeats every 
0.01 sec, which is the period of the blade-passage frequency (100 Hz). The cepstrum 
echo spike at X = 0.002 sec corresponds to a » 0.18, and the echo spike at 
t = 0.020 sec corresponds to a * 0.13; these values of a are considerably differ- 
ent from the expected echo strengths (i.e., 0.49 and 0.10). This discrepancy in the 
data is partly due to the presence of background noise, which we have seen can lower 
the levels of the cepstrum echo peaks. It is also partly because the reflections are 
not perfect duplicates of the direct signals, another effect which tends to lower and 
broaden the echo spikes . 

Some engineering judgment is required to remove the echo spikes since the 
cepstrum levels due to the direct signal at these values of x are not known. The 
edited cepstrum is shown in figure 16(c), and the reconstructed spectrum is shown in 
figure 16(d). The harmonic ripple is almost completely removed from the measured 
spectrum except at the higher frequencies. The difference between the measured 
(contaminated) spectrum and the corrected spectrum is shown in figure 16(e). Al- 
though the difference in overall levels is 1 .0 dB, the correction to individual 
frequency levels oscillates approximately +4.0 dB. 

Another example of data for a typical rotor-test condition (i.e., 55 knots, 6° 
rotor angle of attack, and microphone 2, f Bp = 100 Hz) is given in figure 17 for 
another inflow location. The ripple in the spectrum has a period of about 400 Hz 
(for an echo time delay of 0.0025 sec). From the impulsive data results (table III 
of ref. 20), we expect reflections at i = 0.0031 and 0.0098 sec. A large spike is 
shown in the cepstrum (fig. 17(b)) at x = 0.0026 sec. A spike at x = 0.0090 sec 
is also shown, but its amplitude is much less than the one at x = 0.0026 sec. This 
spike at x = 0.0090 sec could also be part of the natural periodicity of the 
cepstrum. If we had no prior knowledge of the possible echoes, this second reflec- 
tion could easily be overlooked in the cepstrum. Thus, some knowledge of the proba- 
ble echo delay time is desirable for the success of this technique. Note again the 
somewhat periodic character of the cepstrum. The effect of removing the first echo 
is shown in figures 17(c) and 17(d), which show the edited cepstrum and the recon- 
structed spectrum. Again, the ripple at the higher frequencies has been virtually 
removed. Although the correction to the overall sound level is 0.1 dB, the spectral 
differences are ±0.6 dB. (See fig. 17(e).) 

Coefficient-averaged spectra .- To investigate the deterministic content of the 
data, a second method of spectral averaging was used. In this case, the Fourier 
coefficients of 50 spectra rather than the power of 50 spectra, are averaged as 
follows: 
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This coefficient-averaging approach serves to average out random contributions to the 
signal, such as the helicopter-rotor broadband noise r(t) and the tunnel background 
noise contamination n(t). The same data of figure 16 were averaged in this manner 
and are presented in figure 18. For the comparison of the two methods of averaging, 
the overall sound level from the power-averaging method is 112.6 dB, whereas the 
overall sound level from the coefficient-averaging method is 107.2 dB. The blade- 
passage frequency, at approximately 100 Hz, has the same sound level for either aver- 
aging method. This indicates that it is primarily a pure harmonic signal with little 
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random content. (Recall that this was already suspected based on the 0.01 -sec peri- 
odicity in the cepstrum of fig. 16(b).) Between 500 and 2500 Hz, the power -aver aged 
levels are consistently larger than the coefficient-averaged levels, typically by 
5 to 10 dB. Above 2500 Hz, the power -aver aged levels are consistently 1 5 to 20 dB 
larger than the coefficient-averaged levels. This indicates that except at the fun- 
damental frequency, the rotor signal contains a larger contribution from random 
sources than from deterministic sources or phenomena, although the character of the 
signal is quite discrete or tonal. 


Discussion 

A comparison of the cepstra obtained from the two averaging methods is interest- 
ing. (See figs. 16(b) and 18(b).) The cepstrum from the coefficient-averaged spec- 
trum (fig. 18(b)) is much "noisier" than the power-averaged cepstrum (fig. 16(b)) 
and displays many spikes. This is probably due to the more tonal nature of the 
coefficient-averaged spectrum. A prominent spike is again shown at x = 0.002 sec 
and a less prominent spike is shown at t = 0.020. Interestingly, the cepstrum peak 
amplitudes at these echo times are almost the same in the two cepstra, and this 
indicates that the two averaging methods retain the same echo characteristics. 

The fact that the cepstrum echo spikes have the same sound level for both aver- 
aging methods indicates that the noise N(f) present in the power-averaged spectrum 
does not have a measurable effect on the total cepstrum at the echo peak location. 

It appears that the factor which has lowered the cepstrum echo peaks from the ex- 
pected sound levels is the nonlinearity of the reflection process. The conclusion 
that echo distortion effects are more important than additive noise effects was also 
reached from the use of the power cepstrum with hydroacoustic data containing echoes 
(ref. 3) and in a general study of the comparative effects of distortion and noise on 
the cepstrum (ref. 22). Ideally, the removal of any nonlinearity in the echo trans- 
fer function is desirable, as illustrated by the simple example in section 4. Addi- 
tional work was performed in an attempt to estimate a transfer function to correct 
the measured rotor data; however, this work proved largely unsuccessful. The lack 
of success is attributed to two factors: (1) the rotor is a large (i.e., 2.8-m- 

diameter) distributed-noise source which could receive echoes from reflectors dif- 
ferent from those used to calculate the applied H(f), and (2) an accurate H(f) is 
difficult to calculate from the amount and type of available data from this experi- 
ment. In spite of these problems, it is still concluded that editing of the power 
cepstrum provides an improved estimate of the rotor free-field spectrum. 


7. CONCLUSIONS 

The application of the power cepstrum to remove acoustic reflections from mea- 
sured helicopter-rotor acoustic spectra has been investigated. It was found that a 
previously applied correction to the reconstructed spectrum, which would require 
prior knowledge of the echo amplitude or path length, is not required or correct. 
Although information of the echo amplitude and delay time is helpful in applying this 
technique, the technique can often be applied with no prior information of the echo 
characteristics, it was also found that for the case of an exact echoed signal, the 
amplitude of the cepstrum echo spike at the echo delay time is linearly related to 
the relative amplitude of the echo in the time domain. It was found that for the 
best application of the technique, the entire measured spectrum must be from the 
source signal. If part of the source signal is masked by any extraneous non-source- 
related signal, the echo ripple in the power spectrum is not well represented. In 


21 


this case the cepstrum does not yield the desired echo characteristics and aliasing 
in the cepstral domain can occur because of the altered effective sample rate in the 
frequency domain. It was also found that to prevent aliasing of the echo delay time 
in the cepstral domain, the spectral analysis bandwidth (the inverse of the time 
series record) must be less than one-half of the echo ripple frequency (the inverse 
of the echo delay time). 

The power cepstrum was found to be useful for removing some of the contamination 
from acoustic reflections in measured helicopter-rotor acoustic spectra. The cep- 
strum editing process was limited by the nonlinear transfer function acting on the 
acoustic reflections. Although cepstrum editing yielded an improved estimate of the 
free-field spectrum, there still remained some error because of the nonlinear trans- 
fer function of the echo process. An alternate procedure was proposed which would 
allow the complete correction of a contaminated spectrum by use of both the echo 
transfer function and the echo delay time. Although this procedure would remove the 
need for editing in the cepstral domain, it requires a knowledge of the transfer 
function of the echo process. 


NASA Langley Research Center 
Hampton, VA 23665-5225 
March 10, 1986 
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APPENDIX A 


AMPLITUDE OF ECHO CEP STRUM AT T = 0 
The derivation in this appendix is intended to prove that 

ln( 1 + a 2 ) - a Q = 0 


A more complete expression of equation (10), the power cepstrum of g(t) with one 
ideal echo, is helpful to this derivation and is given as follows: 
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_f N “ f N 


D 


2. i2irfx 


/.:■ 

N 


+ 41” cos(2irft ) e l2irfT df 
T ■ e 


- — d> 2 f N — [l + cos(4TTft )] e 

2 J.f 2 e 

N 


i2irf T 


df 


J 4> 3 j' N j[ 3 cos(2uft e ) + cos(6irft e )] e 
-f N 


i2trf x 


df 


_ 1 * 4 f N -l-[3 + 4 cos(4irft ) + cos(8nft )] e 
4 18 e e 

N 


i2*irf T 


df 




( 2irft ) + 5 cos(6irft ) + cos(10irft )] e l2lTfT df 


cos 


- i ik 6 f N — f 1 0 + 15 cos ( 4irf t ) + 6 cos( 8nft ) + cos(12irft )] e 
6 J 32 L e e e 

N 


i2irfT 


df 


i2TTf T 


■1 <t> 7 f N l-[35 cos(2TTft ) + 21 cos( 6irft ) + ...] df 

7 J f 64 L e e 


(Equation continued on the next page) 
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- ^ <t>° J N 7^8 [ 35 + 56 c °s(4TTft e ) + 28 cos ( 8 7rf ) + ...] e 


i27rf x 


1 9 f^N i 

^ <f> J ~ 256 ^ 2 ^ cos ( 2 irf ) + 84 cos^Trft^) + • ..] e 


i27rf x 


1 JO f f N 1 r 
10 * J 51 2 *• 


126 + 210 cos(4irft ) + 1 20 cos(8irft ) + ...1 e 

<=> J 


i2uf t 


+ ... 


This expression may be reduced to 


C (t) = C (x) + 2f [ln(1 + a ) - a 

P pg N L i 


sin( 2irf x! 

1 — 

0 J 2irf t 
N 


sin [ (nt - x ) 2 irf ] sin[(nt + x)2TTf 1 

L e N J L e N J 


+ f Y a < L e AM L e 

N n=1 n l (nt e ' T)2irf N (nt e + x)27rf N 


The form of the coefficient a Q may be shown to be the series of factors multiplying 

all the integrals of e l2lTfT in equation (A1 ) . Thus, a Q may be written as 

- _ 1 1 a 2 . 1 2 A . 1 10 .6 . 1 35 ,8 1 1 26 J 0 ^ 

0 22^ 48 + 6 32 ^ + 8 128 ^ + 10 512 * 


1 , 2 , 3 ,4 5 6 35 ,8 63 JO 

4 ^ + 32 < * > + 96 ^ + 1024 *** + 2560 ^ + 


The entire factor multiplying the (sin x)/x term at x = 0 is ln(1 + a 2 ) - a f 
To prove that 


ln(1 + a ) - a Q = 0 
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it is established that ln(1 + a 2 ) is equal to the same infinite series as the coef- 
ficient a q* To show this we define the quantity (1 + a 2 ) as a function of <J> using 
the simplifying expression 




2a 


1 + a 


(A2) 


which may be rewritten as 


<j>a - 2a + <j> = 0 


This is a quadratic equation in a with roots 


„ 2 J 4 - 4* 2 , t , J— 

a “ 2* “ ♦ * ♦ Vl “ 


Now, since 0 < a < 1.0, we must take only the root 


a 





(A3) 


From this relation and 


(l _ \fi - <() 2 )( 1 + /1 - <f> 2 ) = 


equation (A3) may be expressed in terms of $ as 


a 


1 + 



( A4) 


To find the quantity 1 4* a 2 in terms of <j), equation (A2) is inverted and 
equation (A4) is used to define a: 


1 


+ a 


2a 

♦ 
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Taking the natural logarithm and separating terms results in 


1 + yi " ♦ 


= In 2 - ln\ 


From reference 26 , we find the useful identity 


(< ^ 0 v , x k (2k - 1 )! 2k , 2 ^ „ 

i\1 + VI + x / = In 2 - 2-r (“1 ) — zr- x (x < 1 


2 2k dcl) 2 


Letting x = i<(> to introduce a negative sign under the square root, we may write the 
identity as 


,(i Wi - f 2 ) - m 2 - f <-n k <M» 2k 


2 2k w> 2 


ln 2 - r (2k - '» 

k=1 2 21c ( k ! ) 2 


By substitution of equation (A6) into (A 5 ), 


, ,, 2. , „ r „ v> <2k - 1 )! .2k 

ln(1 + a ) = ln 2 - In 2 - 2- — 7- <(> 

k=1 2 (k!r 


E ( 2k - 1)1 2k 

"2k 2" * 

k-1 2^ K (k ! ) Z 


Since 


E ( 2k - 1 ) ! ,2k 

"2k 2 * 

k=1 2 (k!) 


1 ,2 3 ,4 5 ,6 35 x 8 63 10 , 

4 * + J 2 * + 96 * + T 024 <,) + ^ 56 ?* + 


this shows that 


ln(1 + o ) - a Q = 0 


Thus, there is no (sin x)/x term because of the echo in the power cepstrum at 

T = 0. 
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APPENDIX B 


AMPLITUDE OF ECHO CEPSTRUM AT T = 

The derivation in this appendix is intended to prove equation (15): 


a 1 * 



\ (2j - 1)1 

i; 2 2j - 2 (j - Dljl 


♦ 


2j-1 


2a 


In this series, the first fraction (1/(2j - 1)) results from the logarithm expansion 

and the second fraction ((2j - 1)!/[2 2 ^“ 2 (j - 1)!j!]) results from the expansion for 
powers of the cosine function, as discussed in section 2. 

For values of <j> < 1 .0, the series for a 1 converges to the quantity 2a. To 
prove this, we start by rewriting 2a in terms of using equation (A3): 


2a = fG " t - 4> 2 ) 


( B1 ) 


The problem now is to prove that equation (B1 ) converges to equation (15): 


2 oo 2*i 

11 V ( 1 ) (2j - 2)1 

♦ j=i ' 2 ' <3 - 1)! ^ ! 


|(i - 


( B2 ) 


The right side of equation (B2) may be rewritten as a series, with the binomial 
expansion for - x 2 , as 


(, . it ~T5) = 


2 4 6 


2 4 6 8 


<t> 8 + 


(B3 ) 


With the following definitions. 


(- 1 )!! = 1 
( 2 j + 1 ) ! ! = 1 • 
( 2 j ) !! = 2 • 4 • 


3 • 5 ... ( 2 j + 1) 
6 ... ( 2 j ) 
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equation (B3) may be rewritten in the condensed form as 


g . . * 2 ) . t 

j=i 


( 2j - 3)!! A 2j 
— ttt; < P 


(B4) 


By induction it 


can be shown that 


00 

2 £ 


j=1 



( 2 j - 2)1 
(j - 1)!n! 



( 2 j - 3)11 
( 2 j ) i! 


♦ 


2 j 


or that 


( 2j - 2)! 


2 j — 1 . 

2 (3 - 


1)!j! 


= ( 2 j - 3)1! 
( 2 j ) ! ! 


(B5) 


as follows: Step 1 is to show that equation (B5) is true for j = 1: 


[ 2 ( 1 ) - 21 1 

2 2(1| -’(, - on. 


[2(1) - 3]!! 
[ 2 ( 1 )]!! 


which results in 


1_ 

2 


2 


Step 2 is to assume equation (B5) is true for j = k and to show it is true for 
j * k + 1 • 


[20c + 1) - 2]! 

2 2lkt,) -W + ,)! 


[2(k + 1 ) - 3] ! 1 
[20c + 1)]!! 


which reduces to 


(2k)! = (2k - 1 ) ( 2k - 3)1! 

2 2k+1 k! (k + 1 )! (2k +2) (2k)!! 
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With equation (B5), the above expression is rewritten as 


200! = (2k - 1 ) ( 2k - 

2 2k+1 k!(k +1)! (2k + 2)2 2k-1 (k 


which readily reduces to 


Thus, 


2 2k (2k + 2) 


2 2k (2k + 2) 


H - t - |(i - V 


j=1 2 2j ‘ 2 (j - 1 ) ! j I 


2 )! 

- 1 ) !k! 


f 2) 

1 - < 1 > / = 


2a 
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SYMBOLS 


n 


a( nf B p) 


C' 

p 


"pg 


BP 


N 

G(f) 

g(t) 

H(f) 

h(t) 


k 

N(f ) 
n(t) 
R(f ) 

*x ( t) 
r ( t) 

T 

t 

t 


coefficients of cepstrum echo peaks for one-echo case (n = 0, 1 , 2, ...) 

amplitude envelope of helicopter-rotor harmonic tones (nf Bp ) 

power cepstrum, inverse Fourier transform of logarithm of power spectrum 

edited power cepstrum 

power cepstrum of g(t) 

frequency, Hz 

blade-passage frequency of helicopter-rotor acoustic signal, equal to 
number of rotor blades times shaft rotation frequency, Hz 

fundamental tone frequency, Hz 

Nyquist frequency, Hz 

Fourier transform of direct signal g(t) 
general function of direct signal 

Fourier transform of transfer function of echo process 

impulse response function of echo process 

imaginary part of complex conjugate 

power of example exponential functions 

Fourier transform of noise function n(t) 

noise function 

ratio of X Q (f) to X d (f) 

autocorrelation function 

helicopter-rotor broadband noise signal 

period of time series used in Fourier analysis, sec 

time, sec 

echo delay time of one-echo signal, sec 
period of tone frequency, sec 

echo delay time of first echo in two-echo signal, sec 
echo delay time of second echo in two-echo signal, sec 
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X(f) 


Fourier transform of x(t) 


X d (f ) Fourier transform of direct signal x(t) 

X e (f) Fourier transform of echoed signal x(t) 

x(t) general function of time 

a amplitude ratio of first echo to direct signal 

3 amplitude ratio of second echo to direct signal 

T quefrency, sec 



Subscripts: 

avg averaged quantity 

r real part of complex variable 

A superscript asterisk * indicates a complex conjugate. 
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(a) Time history. 



(b) Power spectrum. 

Figure 3.- Exponential signal with one ideal echo (o = 0.! 

and t_ = 0.0018 sec). 
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(c) Power cepstrum. 
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(d) Reconstructed power spectrum. 


Figure 3.- Continued 
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(e) Spectral correction 
Figure 3.- Concluded. 
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Figure 4.- Exponential signal with two ideal echoes (o = 0.! 
t 1 = 0.0018 sec, 8 = 0.2, and t, = 0.0068 sec). 
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(b) Power spectrum. 

Figure 5,- Exponential signal with one ideal echo (a = 0.5 
and t = 0.0018 sec) and 60 dB of additive noise. 
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(a) Time history. 
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(b) Power spectrum. 


Figure 6.- Exponential signal with one ideal echo (a = 0.5 
and t_ = 0.0018 sec) and 48 dB of additive noise. 
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(d) Reconstructed power spectrum. 


Figure 6.- Concluded 
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(a) Time history. 
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(b) Power spectrum. 

Figure 7.- Exponential signal for k = 20 000 with one 
ideal echo (a = 0.5 and t p = 0.0018 sec) and 48 dB 
of additive noise. 



Power 

Cepstrum/2f N 



0.000 0.050 0.100 0.150 0.200 

0.025 0.075 0.125 0.175 0.225x 10' 


Quefrency, seconds 


(c) Power ceps tr urn. 



(d) Reconstructed power spectrum. 


Figure 7.- Concluded 



(a) Time history. 



(b) Power spectrum of direct signal. 
Figure 8.- Impulsive signal with one echo. 
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(c) Power spectrum of echoed signal. 
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Figure 8.- Concluded. 
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Figure 9.- Correction of contaminated data of figure 8(a). 
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(c) Reconstructed power spectrum found from correction for 

H(f) and t e . 



(d) Reconstructed power spectrum found from editing 
of power cepstrunw 


Figure 9,- Concluded 



(b) Power cepstrum. 

Figure 10.- Periodic signal without echo. Off-harmonic 
level is source related. fr = 5 Hz. 
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(a) Power spectrum. 



Quefrency, seconds 
(b) Power cepstrum. 

Figure 11.- Periodic signal with one ideal echo (a = 0.5 and 
t e = 0.7 sec) for source-related off-harmonic level. 
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(c) Edited power cepstrum. 
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(d) Reconstructed power spectrum. 


Figure 11.- Concluded 
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(a) Power spectrum (power-averaged spectrum). 
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Figure 16.- Model-scale helicopter-rotor acoustic data for velocity 
of 60 knots, rotor angle of attack of 2°, and rotor thrust 
coefficient of 0.007. f Bp = 100 Hz; microphone 3. 
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Figure 16.- Continued. 
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(e) Spectral correction 


Figure 16.- Concluded 
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(b) Power cepstrum. 

Figure 17.- Model-scale helicopter-rotor acoustic data for velocity 
of 55 knots , rotor angle of attack of 6°, and rotor thrust 
coefficient of 0.007. f Bp = 100 Hz; microphone 2. 
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(c) Edited power cepstrum. 
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(d) Reconstructed power spectrum. 


Figure 17.- Continued 
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Figure 18.- Coefficient-averaged helicopter-rotor data 

of figure 16. 
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